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ABSTRACT 

We have carried out high-precision timing observations of the eclipsing binary PSR 
J2051— 0827 in the 3.3 years since its discovery. These data indicate that the orbital 
period is decreasing at a rate of Pb = ( — 1 1 ±1 ) x 1CT 12 . If secular, this orbital period 
derivative implies a decay time for the orbit of only 25 Myr which is much shorter than the 
expected timescale for ablation of the companion. We have also measured the proper 
motion of the pulsar to be 5±3masyr _1 . Assuming the pulsar is at the dispersion- 
measure distance this implies a very slow transverse velocity t>t=(30±20) kms^ 1 . This 
combination of low velocity and short orbital period argue against formation of the 
system in the standard manner and we discuss the implications for its evolutionary 
history. 
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1. Introduction 

The discovery of the millisecond pulsar PSR J2051- 
0827 (Stappers et al. 1996) as part of the Parkes 
survey of the southern sky for low-luminosity and 
millisecond pulsars (MSPs) (Manchester et al. 1996, 
Lyne et al. 1998) brought to five the number of 
pulsars which are eclipsed by winds from low-mass 
(< 0.1 M0) companions. It is only the second such 
system that is not a member of a globular cluster 
and has the second-shortest orbital period of any 
binary-radio pulsar, Ph=2.38h. The duration of the 
eclipse, ~10% of the orbital period, and the de- 
lays at the eclipse egress and ingress indicate that 
there is material, probably originating from a wind 
driven from the companion, well beyond the compan- 
ion's Roche lobe. Thus, this system joins both PSR 
B1957+20 (Fruchter, Stinebring, & Taylor 1988) and 
PSR B1744-24A (Lyne et al. 1990) as neutron stars 
which arc in the process of ablating their compan- 
ions and possibly becoming isolated MSPs (Ruder- 
man, Shaham, & Tavani 1989, Bhattacharya & van 
den Hcuvcl 1991). 

Pulse times of arrival are extremely sensitive to 
changes in a pulsar's local environment and along 
the propagation path. Timing measurements of PSR 
J2051— 0827 have already shown that there is struc- 
ture in the electron column density in the eclipse re- 
gion (Stappers et al. 1996) and can provide valuable 
information on the nature of the eclipses (Stappers et 
al. 1998a). The eclipsing behaviour and circularity of 
the orbit suggest that relativistic effects are not going 
to have an observable effect on the timing. However, 
measurements of the changes in the orbital period for 
both PSR B1957+20 (Ryba & Taylor 1991, Arzouma- 
nian, Fruchter, & Taylor 1994) and PSR B1744-24A 
(Nice & Thorsett 1996) suggest that the effects of ei- 
ther the tidal influence of the companion or the abla- 
tion process itself may be evident in the timing data. 
If the pulsar has a reasonable transverse velocity then, 
given its distance, ~1.3kpc, the proper motion would 
also be apparent. 

2. Pulsar Timing Analysis 

Since its discovery, observations of PSR J2051— 0827 
have been made as part of our regular high-precision 
timing program at Parkes (e.g. Bell et al. 1997). 
Timing data have been obtained approximately ev- 
ery three weeks, although there are some large gaps 
resulting from unavailability of the telescope. Ap- 



proximately 2300 valid pulse times of arrival (TO As), 
at centre frequencies near 436, 660, 1520, 1940 and 
2320 MHz, have been obtained using a filtcrbank sys- 
tem and/or the Caltech Fast Pulsar Timing Machine 
(FPTM). The FPTM is based on an autocorrclator 
capable of very fast sampling (Navarro 1994). Two 
filtcrbank systems are used; at the two lower frequen- 
cies, a 2x256x0. 125MHz filtcrbank detects orthogo- 
nal linear polarizations, while at the higher frequen- 
cies a 2 x 64 x 5 MHz filtcrbank detects orthogonal cir- 
cular polarizations. Details of the filterbank-based 
observing systems can be found in Manchester et al. 
(1996), and the correlator-based system is described 
in Navarro (1994). 

The filtcrbank data were reduced offline. A mean 
pulse profile for each frequency channel was produced 
by folding at the topocentric period of the pulsar. 
To compensate for the dispersive delay these pro- 
files were transformed to the Fourier domain, phase 
shifted, transformed back to the time domain and 
summed. In the FPTM the down-converted signal is 
2-bit digitised and autocorrelation functions (ACFs) 
computed. Bandwidths were typically 128 MHz and 
32 MHz for observations above and below 1GHz, re- 
spectively. The ACFs were hardware-integrated at 
the apparent pulsar period for either 60 or 90s. Af- 
ter being transferred to a Sun Sparc-20 computer the 
data were calibrated, dedispersed and summed. 

The typical integration time required to detect the 
pulsar is 60-90 s, although scintillation caused by the 
interstellar medium means that much longer integra- 
tion times are sometimes required to obtain an ac- 
ceptable signal-to-noise ratio. Each integration was 
optimized to maximize orbital phase coverage with 
individual TO As, while preserving a minimum signal- 
to-noise ratio of ~12. The observed profiles were then 
cross-correlated with a standard pulse template at 
each frequency to determine the TOA. Separate tem- 
plates were used for the higher resolution correlator 
data and no offsets between the two data sets were re- 
quired to the limit of the best fitting residuals. These 
data were then fitted using the TEMPO timing analy- 
sis package (Taylor & Weisberg 1989) with the JPL 
DE200 ephemeris (Standish 1982). Binary parame- 
ters were fitted using the timing model of Blandford 
& Tcukolsky (1976). 
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3. Parameters 

A binary timing model which included the pul- 
sar position, dispersion measure, rotational period 
and period derivative, orbital period, projected semi- 
major axis of the orbit, and orbital period deriva- 
tive was used to generate the best fitting parameters 
shown in Table 1. Data which lie between orbital 
phases 0.2 and 0.35 were discarded from all fits be- 
cause of the influence of the excess column density 
in the eclipsing region on these TOAs. This range is 
narrower than that rejected by Stappcrs et al. (1996) 
as the eclipse boundaries are now better defined. In 
all fits, both the eccentricity and the longitude of pe- 
riastron were held fixed at zero to determine the re- 
maining parameters. The eccentricity limit was de- 
termined by fitting simultaneously for the eccentric- 
ity and the remaining parameters for various input 
values of the longitude of periastron. 

The residuals for this best fit are shown in Figure 1 . 
There is unmodelled timing noise on timescales rang- 
ing from less than an orbit to a few orbits. The intra- 
orbit variations can be attributed to the low signal-to- 
noise of some profiles. The errors are from the formal 
fit and assumed that the errors were completely ran- 
dom. In practice narrow-band interference introduces 
ripples into the pulse profile which lead to underes- 
timates of the true error. There are no trends ob- 
served in these changes and a sufficiently large num- 
ber of orbits have been observed for orbital character- 
istics to remain unaffected. Changes on the timescale 
of a few orbits may be due to changes in the elec- 
tron column density. Variations on longer time scales 
similar to those seen for PSR B1957+20 (Ryba & 
Taylor 1991, Arzoumanian, Fruchter, & Taylor 1994) 
and PSR B1744-24A (Nice & Thorsett 1996) are 
not present. A formal fit for a period second deriva- 
tive to the timing data of PSR J2051-0827 yields 
P=(4±2) x 10 _29 s _1 , suggesting that mass motions in 
the PSR J2051-0827 system do not greatly affect the 
rotation of the pulsar. 

If the transverse velocity of a pulsar is sufficient, 
then the contribution of the apparent acceleration 
along the line of sight to the derivatives of the or- 
bital and rotational period of the pulsar is significant 
(Shklovskii 1970). This is especially the case for MSPs 
which have rotational spin-down rates smaller than 
those of the long-period pulsars. Corruption of the 
rotational period derivative means all properties of 
MSPs which are derived from it, such as the char- 



acteristic age, t c = P/2P, and the pulsar spin-down 
energy, E p = 4n 2 IP/P 3 (I = 10 45 gcm~ 2 ), are only 
provisional until the velocity, or an upper limit, can be 
determined. Modelling the lightcurve of the compan- 
ion to PSR J2051-0827 has shown that at least 30% 
of the pulsar's spin-down energy may be required to 
heat the companion star (Stappers et al. 1998b); this 
fraction would increase to nearly 50% for a proper 
motion of lemasyr" 1 (v t ^lOOkms -1 ). However, 
as shown in Table 1, our current measurement of 
the composite proper motion of PSR J2051— 0827 is 
5±3masyr _1 . Assuming the dispersion- measure dis- 
tance to the system, this implies a remarkably slow 
transverse velocity of only v t = (30±20) kms -1 . Thus 
the contribution to the period derivative is 3.4 x 10~ 22 , 
or just 3% of the measured period derivative. 

A significant result from these timing observations 
is that the orbital period of PSR J2051-0827 is de- 
creasing at a rate A = (-ll±l)xl0~ 12 . This Pb is 
some two orders of magnitude greater than the con- 
tribution expected from general relativistic effects, 
Pb = (— 3±1) x 10~ 14 , and the possible influence of 
the Shklovskii term is negligible. If the orbital pe- 
riod derivative were constant then the orbital decay 
time would be only 25 Myr. However, observations of 
PSR B1957+20 indicate that its orbital period deriva- 
tive varies, and even changes sign, on quite short 
time scales (Arzoumanian et al. 1994). The best 
fit obtained for a constant orbital period has an r.m.s 
~35/^s, somewhat greater than for the fit in Tabic 1. 

If the orbital period is decreasing, the pulsar should 
arrive at the ascending node earlier than we would 
predict from the constant orbital period model. Data 
presented in Figure 1, except those obtained in 1997 
July, were split into three, approximately equal, time- 
interval groups. After re-assigning the epoch of as- 
cending node to the middle of each data set, they were 
fitted for orbital and spin parameters. The position 
was held constant for all fits as none of the data sets 
spanned a full year, and the proper motion is small. 
The measured phase shifts are found to be consistent 
with those predicted by the measured orbital period 
derivative. Thus confirming the orbital period of the 
system is presently decreasing. There is insufficient 
data to fit for higher-order variations in the orbital 
period. 
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4. Evolution: Past 

The orbital period of PSR J2051-0827 puts it just 
above the upper edge of the low-mass x-ray binary 
(LMXB) period gap. It is one of only four binary 
pulsars with a period less than 4h. It is interest- 
ing to note that of these four binary pulsars, only 
J1910+0004 (Dcich ct al. 1993) shows no strong ev- 
idence of eclipsing behaviour, and it is probably too 
weak to preclude the possibility that it is also an 
eclipsing system. Thus at least three out of the four 
systems are eclipsed, more than would be expected, 
naively, if eclipses in these systems required inclina- 
tion angles close to 90°. This is further evidence that 
there must be material which extends well beyond the 
companion's Roche lobe. The eclipsing phenomenon 
is perhaps inevitable for such short-period binaries. 

The proper motion and dispersion-measure-derived 
distance of PSR J2051— 0827 indicate its transverse 
velocity is very low (30 kms -1 ). Unless there is an 
unusually large radial velocity component, the space 
velocity will also be small. There are presently ten 
MSPs with estimates of their proper motions and they 
indicate a mean space velocity of ^100 km s -1 . This 
is significantly less than the mean velocity of the long- 
period pulsars (e.g. Lyne & Lorimer 1994). Simula- 
tions by Cordes & Chernoff (1997) and Ramachan- 
dran & Bhattacharya (1997) show that such a low 
mean velocity is expected for the MSP population as 
progenitor systems that receive high kick velocities 
will disrupt during the supernova. Thus, the low ve- 
locity of PSR J2051— 0827 is consistent with its prox- 
imity to the disk (z=650pc) but is problematic when 
considering its past evolution. 

In standard evolutionary models, low-mass binary 
and isolated MSPs are descended from LMXBs (Alpar 
et al. 1982, Bhattacharya & van den Hcuvel 1991). 
LMXBs typically have orbital periods of a few days 
before the supernova explosion of the primary. The 
simulations of Tauris & Bailes (1996) show it is al- 
most impossible for such initially compact systems to 
generate binary MSPs with orbital periods less than 
a day and spatial velocities less than 50kms _1 . This 
is because the shell of the exploding star removes a 
large amount of momentum from the system when 
the neutron star is formed, and the binary recoils. 
The closer the system at the time of the explosion, 
the greater the recoil velocity. Thus we expect the 
shortest orbital-period MSPs to have the greatest ve- 
locities. While the transverse velocity, w t ss220kms _1 



of PSR B1957+20 (Arzoumanian et al. 1994), is con- 
sistent with the standard formation mechanism and 
the simulations of Tauris & Bailes (1996), the low ve- 
locity of PSR J2051-0827 indicates that it probably 
did not form this way. 

Any alternative formation process must explain the 
low velocity, short orbital period, and 4.5 ms pulse pe- 
riod of PSR J2051-0827. The low velocity suggests 
the orbital separation at the time of the explosion was 
wide. The only reasonable option requires a massive 
spiral-in since the formation of the neutron star. Mass 
transfer from a more massive to a less massive body 
leads to spiral-in, but it must occur slowly enough 
for the neutron star to be spun up to a millisecond 
period. Thus the companion cannot have been too 
massive. Our ideal progenitor system is therefore an 
intermediate-mass binary with an orbital period of a 
few months and a secondary in the mass range 2-4 
M . There is enough mass for the system to avoid 
disruption when the primary explodes, a sufficiently 
long orbital period to ensure a low runaway veloc- 
ity, the correct mass-ratio for eventual spiral-in of the 
neutron star into the secondary core, and enough time 
to spin up the pulsar. 

This alternative evolutionary process may explain 
the low velocity of PSR J2051-0827 but nominally its 
end product is a MSP and a massive white dwarf in 
a circular orbit. If the remaining core mass after the 
common-envelope phase was sufficiently small, then 
mass transfer to the neutron star may not have be- 
come unstable. The final spiral-in phase would have 
been avoided and the system may have proceeded in- 
stead to a compact LMXB-like phase. This scenario 
is attractive as the evolution of the system would 
then be similar to that discussed by Ruderman ct al. 
(1989) and would naturally explain the present very 
low-mass companion. 

Alternatively, if the neutron star was formed via 
accretion-induced collapse of a massive white dwarf 
(e.g. Bailyn & Grindlay 1990) then the velocity of the 
system would also be lower as less mass is lost dur- 
ing formation of the neutron star. However, it would 
also require the explosion to be symmetric, otherwise 
the system would still be expected to have a large ve- 
locity (Tauris & Bailes 1996) and would require the 
resultant pulsar to ablate most of the mass from the 
companion. 
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5. Evolution: Future 

Initial ablation of the companion in eclipsing and 
isolated systems is thought to occur during an LMXB 
phase when a 7-ray flux may be generated by the in- 
teraction of the neutron star's magnetic field and the 
accretion disk. This flux causes a wind to be driven 
from the companion star by heating its outer layers 
(Ruderman et al. 1989, Kluzniak et al. 1988). Evap- 
oration of the companion may eventually cause the 
accretion to cease and the spun-up pulsar would then 
be able to turn on. The pulsar may then continue 
to ablate the companion through its relativistic wind 
(Ruderman, Shaham, & Tavani 1989). Following the 
simple calculation of Bhattacharya & van den Heuvel 
(1991), we relate the rotational energy of the pulsar 
to the binding energy of the companion to determine 
if the pulsar can ablate its companion. Assuming the 
ablation process is 1% efficient, i.e., 1% of the inci- 
dent pulsar spin-down energy is available to drive the 
wind, we find PSR J2051-0827 will require a further 
10 9 yr to evaporate its companion. This time scale is 
longer than the 25Myr time scale for orbital decay 
suggesting that tidal effects may be more important 
in completely destroying the companion. If MSPs live 
for 10 9 yr or more this time scale for destruction of 
the companion via orbital decay is inconsistent with 
the observed number of eclipsing MSPs and isolated 
MSPs. 

The value for the orbital period derivative is much 
larger than that predicted by gravitational radia- 
tion losses. Mechanisms whereby angular momen- 
tum is lost from the system through anisotropic mass 
loss from the companion's surface (Banit & Shaham 
1992, Eichler 1992, Brookshaw & Tavani 1993) were 
proposed to explain the orbital period derivative of 
PSR B1957+20. However, they require a large mass- 
loss rate which is hard to reconcile with the low elec- 
tron densities measured in the eclipse region of either 
PSR B1957+20 (e.g. Fruchter et al. 1990) or PSR 
J2051-0827 (Stappers et al. 1996). Moreover they 
are unable to reproduce the rapid change in sign and 
magnitude of the orbital period derivative that has 
been measured for PSR B1957+20. 

These variations in the period derivative have been 
likened to the quasi-cyclic variations seen in other 
close binaries where magnetic fields of the companion 
are important (Arzoumanian et al. 1994). Based on 
this idea, Appelgate & Shaham (1994), developed a 
model where the companion's magnetic activity and 



wind generate a torque which prevents it from co- 
rotating and results in dissipation of tidal energy. Al- 
though we have measured only a secular change in the 
orbital period derivative of PSR J2051— 0827 at this 
stage, given its companion mass is similar to that for 
PSR B1957+20, and the orbital separation is less, we 
might expect it too will have a variable orbital period 
derivative. Clearly if such variation in the orbital pe- 
riod derivative were measured for PSR J2051— 0827, 
it would indicate it will live longer than 25Myr and 
thus alleviate the birthrate problem. 

The low electron column densities measured in the 
eclipse region for PSR J2051— 0827 suggest that the 
current evolutionary phase is a long-lived one. How- 
ever, the large orbital period derivative, if it is secu- 
lar, provides a mechanism to reduce this system to an 
isolated MSP. The low velocity of the pulsar is a pow- 
erful constraint on the progenitor system and favours 
its having evolved from an initially intermediate mass 
binary. 

We thank S. Johnston and the referee for helpful 
comments on the manuscript. We are grateful to S. 
Anderson, M. Britton, A. Hughes, V. Kaspi, A. Lyne 
and the Parkes Observatory staff for assistance with 
observations. B. W. S. received support from an ANU 
PhD scholarship and the ATNF student program. 

REFERENCES 

Alpar, M. A., Cheng, A. F., Ruderman, M. A., & 
Shaham, J. 1982, Nature, 300, 728. 

Applegate, J. H. & Shaham, J. 1994, ApJ, 436, 312. 

Arzoumanian, Z., Fruchter, A. S., & Taylor, J. H. 
1994, ApJ, 426, L85. 

Bailyn, C. D. & Grindlay, J. E. 1990, ApJ, 353, 159. 

Banit, M., Ruderman, M. A., Shaham, J., & Apple- 
gate, J. H. 1993, ApJ, 415, 779. 

Banit, M. & Shaham, J. 1992, ApJ, 388, L19. 

Bell, J. F., Bailes, M., Manchester, R. N., Lyne, A. C, 
Camilo, F., & Sandhu, J. S. 1997, MNRAS, 
286, 463. 

Bhattacharya, D. & van den Heuvel, E. P. J. 1991, 

Phys.Rep. , 203, 1. 
Blandford, R. & Teukolsky, S. A. 1976, ApJ, 205, 

580. 

Brookshaw, L. & Tavani, M. 1993, ApJ, 410, 719. 
Cordes, J. M. & Chernoff, D. F. 1997, ApJ, 482, 971. 
Deich, W. T. S., Middleditch, J., Anderson, S. B., 

Kulkarni, S. R., Prince, T. A., & Wolszczan, 

A. 1993, ApJ, 410, L95. 
Eichler, D. 1992, MNRAS, 254, IIP. 



5 



Fruchter, A. S. et al. 1990, ApJ, 351, 642. 
Fruchter, A. S., Stinebring, D. R., & Taylor, J. H. 

1988, Nature, 333, 237. 
Kluzniak, W., Ruderman, M., Shaham, J., & Tavani, 

M. 1988, Nature, 334, 225. 
Lync, A. G. & Lorimer, D. R. 1994, Nature, 369, 127. 
Lyne, A. G. et al. 1990, Nature, 347, 650. 
Lyne, A. G. et al. 1998, MNRAS, . in press. 
Manchester, R. N. et al. 1996, MNRAS, 279, 1235. 
Navarro, J. 1994. PhD thesis, California Institute of 

Technology. 

Nice, D. J. & Thorsett, S. E. 1996, in Pulsars: Prob- 
lems and Progress, IAU Colloquium No. 160, 
ed. M. Bailes, S. Johnston, & M. Walker, As- 
tronomical Society of the Pacific, p. 523. 

Ramachandran, R. & Bhattacharya, D. 1997, MN- 
RAS, 288, 565. 

Ruderman, M., Shaham, J., & Tavani, M. 1989, ApJ, 
336, 507. 

Ruderman, M., Shaham, J., Tavani, M., & Eichler, 

D. 1989, ApJ, 343, 292. 
Ryba, M. F. & Taylor, J. H. 1991, ApJ, 380, 557. 
Shklovskii, I. S. 1970, Soviet Astron., 13, 562. 
Standish, E. M. 1982, A&A, 114, 297. 
Stappers, B. W., Bailes, B., Lync, A. G., Camilo, F., 

Manchester, R. N., Sandhu, J. S., Toscano, M., 

& Bell, J. F. 1998a, . In prep. 
Stappers, B. W. et al. 1996, ApJ, 465, L119. 
Stappers, B. W., van Kerkwijk, M. H., Lane, B., & 

Kulkarni, S. R. 1998b, . In prep. 
Tauris, T. M. & Bailes, M. 1996, A&A, 315, 432. 
Taylor, J. H. & Weisberg, J. M. 1989, ApJ, 345, 434. 



Table 1: Parameters of PSR J2051-0827 



Right Ascension (J2000) 


20 h 51 m 07!5130(2) 


Declination (J2000) 


-08°27'37'.'782(6) 


Proper motion in R.A. 


1(2) mas yr^ 1 


Proper motion in Deck 


—5(3) mas yr _1 


Epoch of Period (MJD) 


49530.0 


Period (s) 


0.0045086417433540(7) 


Period Derivative (xl0~ 20 ) 


1.272(2) 


Dispersion Measure (cm~ 3 pc) 


20.7458(2) 


Orbital Period (days) 


0.0991102650(2) 


Orbital Period Derivative (xl0~ 12 ) 


-10.8(10) 


Projected Semi-major Axis (lt-s) 


0.045076(1) 


Eccentricity 


< 8 x 10~ 5 


Epoch of Ascending Node (MJD) 


50015.422291 


R.M.S. timing residual (/is) 


32 


Distance (kpc) 


1.3 


Galactic Longitude (degrees) 


39.19 


Galactic Latitude (degrees) 


-30.41 


Mass function (M Q ) 


1.00108(7) xl0~ 5 


Minimum Companion mass (Mq) 


0.027 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v4.0. 
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Fig. 1. — Post-fit timing residuals for ^2300 timing measurements of PSR J2051— 0827 at a number of different 
frequencies (see text) for the parameters in Table 1. 
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